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A model field theory, in which the interaction between quarks is mediated by 
dressed vector boson exchange, is used to analyse the pionic sector of QCD. It is 
shown that this model, which incorporates dynamical chiral symmetry breaking, 
asymptotic freedom and quark confinement, allows one to calculate f n , m n , r n and 
the partial wave amplitudes in tt-tt scattering and obtain good agreement with the 
experimental data, with the latter being well described up to energies E ~ 700 MeV. 
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I. INTRODUCTION 



Chiral Perturbation Theory (ChPT) has been applied extensively and successfully to a broad 
range of low energy strong interaction phenomena; both purely mesonic and with the inclusion 
of nucleons 0. A fundamental ingredient of ChPT, and, indeed, what it emphatically establishes, 
is the notion that the theory underlying the strong interaction must manifest dynamically broken 
chiral symmetry. 

Fundamental to obtaining physical information from ChPT, in the sense of calculating ex- 
perimental observables, is the fact that at first nonleading order [0(-E 4 )]; i.e., to fourth order in 
derivatives of the chiral field, one obtains a complete effective Lagrangian by including one-pion- 
loop graphs which introduce ten counterterms (whose coefficients are infinite). Higher order loops 
do not modify the Lagrangian at this order Q but the ChPT La grangian is only completely defined 
upon the inclusion of these one-loop contributions. 

This framework can be used to analyse the pionic sector of the strong interaction: in particular, 
7r-7r scattering ; and it is quite successful. One can choose the coefficients of the counterterms 
in such a way as to obtain rather good agreement with the partial wave amplitudes. The approach 
of Ref. H augments standard ChPT by the inclusion of scalar and vector resonances and obtains 
good agreement with the experimental S-, P-, D- and F-wave scattering lengths and reasonable 
agreement with the phase shifts up to E ~ 800 MeV. 

There is another approach to modelling QCD which we wish to employ to study the pionic 
sector of QCD, in general, and tt-tt scattering in particular. This approach, as will be seen, 
actually contains ChPT in the sense that it generates an effective Lagrangian for the chiral field 
with the same structure as that employed in ChPT but with the difference that the coefficients of 
each term are finite with values determined by three parameters associated with the quark-gluon 
substructure of the model. Our approach is based on a model field theory, with elementary quark 
fields interacting via dressed vector boson exchange, that has been referred to as the Global Colour- 
symmetry Model (GCM) and which was first described in Ref. Q and extended to the 7r-sector in 
Ref. 0. The extension to other mesons and to the baryon is reviewed in Ref. ||. 

Including quarks explicitly, in the sense of actually carrying out the operations on the generating 
functional for the model field theory, Z, that lead to an expression of Z in terms of the chiral field 
so that the relation between the coefficients in the effective action and those describing the model 
quark-gluon dynamics is made explicit, is one property that sets this approach apart from ChPT. 
Quark confinement is incorporated by ensuring that the quark propagator obtained in the model 
has no singularities that can give rise to absorptive parts associated with free-quark production in 
the Feynman diagrams for any given process. (This is a sufficient but not necessary condition for 
quark confinement, as discussed in Ref. Q.) 

In the GCM the effective action for mesonic degrees of freedom is generated by a derivative 
expansion of the fermion determinant, that arises through an Hubbard-Stratonovich transformation 
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and subsequent integration over the quark fields, in a manner analogous to that employed in the 
Nambu-Jona-Lasinio model (NJL) flOfl . In contrast to the NJL studies, however, each fermion 
loop that arises is finite because the meson fields are not pointlike but, rather, they are extended 
objects whose internal structure is given by Bethe-Salpeter amplitudes obtained via solution of the 
ladder Bethe-Salpeter equation. Jll[] These Bethe-Salpeter amplitudes play the same role as the 
phenomenological form factors employed in the calculations of intermediate energy nuclear physics 
but, it should be stressed, they are calculated quantities in the GCM. 

Another property that sets the GCM apart from ChPT is the fact that meson loops are also 
finite. This is because the meson-meson vertices are not pointlike; i.e., there are no contact 7r(x) 4 
interactions, but instead each interaction has a momentum dependence determined by an under- 
lying quark loop structure. This quark substructure generates "form factors" at the vertices and 
therefore the tree level effective action only receives finite corrections from one- loop diagrams, a 
fact which is illustrated in Ref. [0|. This is an important difference from ChPT. The studies of 



Ref. [ 13 1 also suggest that, in our model, the the loop corrections are small near threshold. 

The GCM is a model which microscopically manifests dynamical chiral symmetry breaking 
(DCSB) via the mechanism well know from phenomenological Schwinger-Dyson equation (SDE) 
studies |l4j and hence it will have in common with ChPT all of the consequences entailed by this. 
Herein we use the phenomenon of tt-tt scattering to illustrate this. Of course, in specifying a model 
one has lost the "model independence" which is a beauty of ChPT. However, one might argue 
that a stronger connection with QCD has been obtained since, in addition to DCSB, the GCM 
incorporates the ultraviolet renormalisation group features of QCD by construction. This enables 
one to connect observable quantities with the dynamical mass scale of QCD, Aqcd ; for example, 
and also to study their dependence on the qualitative and quantitative features of the propagators 
of the elementary excitations in QCD. 

In Sec. II we briefly summarise the GCM and the procedure that allows one to proceed to 
an effective action for the Goldstone modes in the model. The action one obtains is of the same 
structural form as would be obtained in the NJL model but the calculation of physical observables 
is quite different due to the internal structure of the mesons. In this section we also present the 
expressions for /„-, m n and r n in the GCM. In Sec. Ill we study tt-tt scattering at tree level. We 
calculate the scattering amplitude, A(s,t,u), and the five lowest order partial wave amplitudes. In 
addition, through exploratory parameter fitting, we demonstrate that the GCM can reasonably be 
expected to describe this sector of QCD rather well. In Sec. IV we proceed to use the GCM to 
calculate the quantities that characterise the pion sector. We present a three parameter version 
of the model, the parameters characterising the nature of the quark-quark interaction at small k 2 , 
and show that a rather good description of the scattering lengths and partial wave amplitudes is 
easy to obtain. We summarise our results and present our conclusions in Sec. V. 
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II. EFFECTIVE ACTION FOR PIONS 



The action for the GCM is, in Euclidean metric (with {7^,71/} = 25^ and 5^ u the Kronecker 
delta): 



S[q,q} = J AgWlT^ + Mj^ + ^AA^Wg'^^-^fe) ■ 
Here M is the quark mass matrix, j^(x) = q(x)^- r y ll q(x) and 

g 2 D^ lu (x ~y) = J 



f d 4 k ik 




x k^k u 


1 (27r) 4e 


k 2 \ 


^ k 2 



(1) 



(2) 



is the dressed gluon propagator in the model. The specification of Fi(k 2 ) completes the definition 
of the model. (It should be noted that the action for an Abelian theory can be expressed exactly 
in this form, with F^(k 2 ) = 1, after functional integration over the gauge field. fi~5|| ) 

The generating functional for the field theory associated with this model action can be written: 



z bhV\ = J VqDq exp -S[q,q) + / 



S[q,q]+ / d x [rj(x)q(x) + q(x)rj(x)] 



(3) 



where rj and r\ are Grassmannian sources for the fermion fields. Following the bilocal bosonisation 
procedure described in detail in Ref. |@], a crucial step in which is the application of the bilo- 
cal analogue of the Hubbard-Stratonovich transformation familiar from statistical mechanics, this 
generating functional can be rewritten as: 



Z[rj,r]] = J VU exp [— S[U] + source terms] 



which now involves the local chiral field 



U(x) = exp i 



(4) 



(5) 



where {T;i = 1, . . . , iV f 2 - 1} are the analogue of the Pauli matrices in SUf(N\ 
The action in Eq. (Q) is 

S[U] =n[U]+iT[U] = -TrLn [g-\x,y;[U])] + I[U] 

where 



fj- 



(6) 



g-\x, y] [u]) = 

G- 1 (x-y;U) + \P R 



+ Pl 



x + y 



B(x - y) 



G- 1 (x-y;U)=-fdA(x-y) + {p R U + P L W}B(x-y) , 



(7) 
(8) 



and 
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I[U] 



d^w 



1 



75tr 



8N L 



d 4 ztr 



[P R U + P L tf) B{-z)T (tr [r'G(z; U) 



+ 



+ 



1 



d 4 ztr 



-fM-z;U) 8^ A{z)-5\z) 



is the term necessary to complete the square in the Hubbard transformation. In Eq. (|9|) 

= I (I + 75) , Pl = 1 (I 



75 J 



(9) 



(10) 



2 v ,u, , 2 

are the right and left handed helicity projection operators. 

The vacuum or ground state configuration of the chiral field is, as usual, defined as that 
configuration U for which 

4fS[U] . = • (11) 

In the present case Eq. ( |Tl|) entails that G~ 1 (z; U) is the solution of the ladder Schwinger-Dyson 
equation which has the following form in momentum space: 

d 4 k 
(2^ 

If the fermions are massless; i.e., M = 0, then this (vacuum) equation admits the same solution 
for every spacetime independent, unitary flavour matrix U. This signals that chiral symmetry is 
dynamically broken if the ladder SDE admits a nonzero solution for B(p 2 ); a well known result 
seen here from another perspective. 

We analyse the action in Eq. (|j) using a derivative expansion of the fermion determinant 
(TrLn). In this process we expand about the vacuum configuration U and, with the form of I[U] in 
Eq. @) , our "bare" pion fields are not pointlike but actually have an internal structure described by 
the ladder Bethe-Salpeter equation pd}| : the Bethe-Salpeter amplitude and fermion mass function, 
B(p 2 ), are identical because of the identity between the ladder SDE and Bethe-Salpeter equations 



G- 1 {p ] U) = i 1 -p + M + -J —- g 2 D^(p - k)j„G(k; U) lv 



(12) 



in vector exchange theories with dynamically broken chiral symmetry. [12| We have already made 
use of this result in writing Eqs. (0-§). 

The expansion of the TrLn generates terms whose coefficients are constructed from quark loops 
coupling to varying numbers of external pion fields; for example, the diagram in Fig. 1 is associated 
with the expression for f n . In this, and the diagrams associated with the other coefficients, the 



quark lines represent dressed quark propagators obtained from the solution of Eq. (12) and the 
pion Bethe-Salpeter amplitude provides an intrinsic cutoff on the momentum integration. As 
a consequence, each of the coefficients is finite and the action is completely determined, with no 
arbitrary finite or infinite constants being necessary, once the dressed gluon propagator is specified. 
[The only exception to this is the fermion condensate which, in the more sophisticated versions 
of the GCM (and QCD itself), diverges with the upper bound on the momentum integral in a 
manner governed by the anomalous dimension of the fermion propagator so that m q 2 {qq) u '2 is a 
renormalisation point invariant. This will become clear below.] 
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The real part of the action in Eq. (||) involves even numbers of pions and, to 0(E 4 ), it is @: 



n[u] 



<fx I ^tr 



4 



iV c ^itr 



+ 7V c K 2 tr 



(2/ f -U -U ] 



(13) 



N c K 3 -ti IdMd^dMd^ 



where = diag(m^ 2 ) ^^2) is the quark mass matrix (p^M^ is renormalisation point invariant 
in QCD and in the more sophisticated versions of the GCM) and 

U{x)=et f -^ x) . (14) 

At this point we have specialised to the case of SU{(2) [f are the Pauli matrices] and we have chosen 
the vacuum configuration U = if which preserves SU^ (2) symmetry and necessarily follows when 
M 7^ 0. The structural form of this part of the action is the same as that of Ref. flofl , 

The imaginary part of the action involves odd numbers of pions and the leading order term is 
0(£ 5 ): 



iT[U] 



d*x 



Nr. 



1 



240vr2 e ^ 2 tr 



(fPU{x) - U\x)fj\ d tl U{x)d v U(x)d x U{x)dpU{x) 



(15) 



At lowest order in the pion field [0(7r 5 )] this term yields the Wess-Zumino term |0| and when 
vector mesons are included this imaginary part of the action contains all of the "anomalous" 



interactions. [17] 

In the GCM, as mentioned above, the coefficients in Eq. (??) (/„-, p, etc.) are completely deter- 
mined once the fermion propagator is known and this follows as the solution of the fermion SDE, 



Eq. (12), once the model gluon propagator, Eq. (|2|), is specified. Writing the fermion propagator 
in the form 



then: 



G(p) = -il'pavip) + °s{p) 



(16) 



J TT 



^2 J dxB2 (°V ~ 2 [ a so's + xa v a' v ] - x a s cr' s - (a' s ) 



J \ 2 



'vf 



(17) 



(18) 



and Ki, K2, K% are given in the appendix. 

One can also use this approach to study the electromagnetic pion form factor, F n (q 2 ), and this 
will be described elsewhere. [|l8| From the generalised impulse approximation diagram, Fig. (|2|), one 
obtains the piece of the charge radius that is regular in the chiral limit, m n — > 0, and only weakly 
dependent on m^. Preliminary studies suggest that this is the dominant piece at the physical value 



of m n jl9|, a result supported by the lattice QCD studies of Ref. |20|. [In ChPT r n ~ lnm 2 as 
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m n — > 0. Since the pion is a point particle in ChPT it derives its form factor from self dressing; 
the lowest order process being interaction with the photon via a pion loop arising from a tttt — > 7T7t 
interaction. The loop will, of course, have a branch cut at the mass of the loop particle - in this case 
m n - and this is the origin of the logarithmic divergence of r w in ChPT. [21] The preliminary studies 



of Ref. [ 19 1 suggest, however, that at the physical value of m n the charge radius is dominated by 
the regular contribution from Fig. (pi) which is not present in ChPT.] 



III. tt-tt SCATTERING AT TREE LEVEL 



A. Tree level tt-tt scattering amplitude: A(s, t, u) 



In considering tt-tt scattering we need only consider the real part of the action in Eq. (??). It 
will be observed that this differs from that of Eq. (5) in Ref. in SU{(2) the a\ and 02 terms 
therein are contained within the ^1,2,3 terms in our action which actually has a richer structure. 
This is easily seen using the relations: 



tr 



tr 



(duUdprt 



and 



tr 



d^UdvrfdnUdvlfi 



1 

2 tr 



dpUdvUi] tr faudvrf 



(19) 



(20) 



2daSdaSd v v l d u v l - 2d.sd u sd ll v i d u v i + d.v'd.v'd.v'd^ - d^d^d^d^v 



where 



S( x ) = -tr [U{x) + U\x) 



V(x) 



1 

4 tr 



(U(x) - U\x] 



(21) 
(22) 



As usual, the tt-tt scattering kernel can be written as follows (with s, t, and u the usual 
Mandelstam variables): 



T a /3;-yS = A(s, t, u) SapS^s + A{t, s, u) 5 ai 5p S + A(u, t, s) 5 aS 5^ 

which can be decomposed into amplitudes with definite isospin as follows: 

T°(s, t,u) = 2 T + _ ;+ _ - 2 T +0;+0 + T 00; oo = 3A{s, t, u) + A(t, s, u) + A(u, t, s) 

^(s^u) = 2T +0;+0 - T ++;++ = A(t,s,u) - A(u,t,s) 

T 2 (s, t,u) = 2 T 00 ;oo + 2 T 00 ;+- + 2+ +;++ = A(t, s, u) + A(u, t, s) . 



(23) 



(24) 
(25) 
(26) 



The partial wave amplitudes are obtained from these expressions as follows: 
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TRs) = J- J^dyP^T^s^u) . (27) 

(When the pions are on mass shell one has, of course, s + t + u = 4m 2 and, in evaluating the partial 
wave amplitudes it is simplest to work in the centre of mass frame where t = — (s — 4m 2 )(1 — y) /2.) 
To calculate A(s,t,u) it is necessary to expand the action of Eq. (??) to 0(tt(x) ): 



f d 4 x < 






.2 



8^%^ + m 2 7rV 



lllf (vrV) 2 
24 /2 1 ' 



f 1 

J IT 



where 



which is renormalisation point invariant. 

In order to facilitate comparisons with Ref. M we have defined 

ai = iV c (2^ a - 2K 2 - K 3 ) , 
a 2 = 2N C K 3 

in Eq. (|H). 

It is now a simple matter to proceed from Eq. (|2"8|) and obtain 



.A(s, i, u) 



mi + 2 s — t — u 



3/ 2 

Ki (-12ml + 6m 2 n (s + t + u) + 2 s 2 - t 



2 u 2 - 2(st + su + 



4iV c 

+ivT 2 (-2m 2 +s) (-2m 2 +i + u) + K 3 (-2ml + m 2 n (s + t + u) - tu 



(28) 



(29) 



(30) 
(31) 



(32) 



[It will be remembered that Eq. (^) is a Euclidean space action. In deriving Eq. ( |32| ) we calculated 
an algebraic expression for A(s^,t^,u^) and used this to define an analytic continuation to the 



physical (s,t,u) domain.] It is clear that the Weinberg result p3] is contained within Eq. (32): one 
must simply set K \ = = K 2 = K 3 ; i.e., ignore the terms at 0(E 4 ). 



B. Tree level tt-tt partial wave amplitudes 

The expressions in Eq. (|26|) can now be used to obtain the scattering amplitudes at tree level 
and here we present our calculated expressions for the first five. All of these expressions have no 
imaginary part because we are working at tree level: 



S 



I 



1 



7ml + 



8m 



4 r 



327T/2 V * ft 

Am 2 

J IT 



32 

5(ai + a 2 ) + y^c^i 



+2x 



1 + ^ [4ai + 3a 2 + AN c Kx] + [H<*i + 7a 2 ] 



96 7T/2 
-1 

32V/f 



1 + 



2m" 



4m 



2 r 



(2 
J - 



(a 2 - 2ai + 4A^ CJ PCi) + — (a 2 - 2ai) 

/ 7T 



(33) 
(34) 



8m; 



ai + a 2 - 2 N cKi 



+x 



1 



4m 



f 2 

(a j + 2a 2 ) x 2 
240tt/4 : 

(2a^ + ag) x 2 



£ [2ai + 3a 2 + 4^^] 



3/< 



[ai + 2a 2 ] 



(35) 
(36) 
(37) 



4807T/4 ' 

where x = s — Am 2 and hence threshold corresponds to x = 0. The violation of unitarity by these 
tree level amplitudes will become important as one moves away from threshold. 
Near threshold the partial wave amplitudes have the form: 



Re 



T l I (u)]=u l (aj + ubj + Oiu 2 )), 



(38) 



where u = x/(2m 7T ) 2 , and a first test of the GCM in this sector is the evaluation of these threshold 
parameters. The dimensionless scattering lengths 



-01 



lu=0 ' a l — du 1 ! 



u=0 



! "0 



T, 



\ U =Q ' fl 2 _ S7 7J -2 







d 2 jiO 



u=0 



a\ 



d 2 rp2 

du 2 2 



u=0 



(39) 



are not influenced by the imaginary part of the partial wave amplitudes and hence they provide a 
useful, simple first test of tree level calculations in a given model. 

The best known of the scattering lengths are a§ and a 2 , and in Ref. it is argued that, 
experimentally, 



0.20 ±0.01, al = -0.037 ± 0.004. 



(40) 



The others are rather less well known. In Ref. [25], however, based on the experiments of 



Ref. [26 1, the following values are presented: 

a° = 0.26 ± 0.005, a 2 , = -0.028 ± 0.012, a{ = 0.038 ± 0.002, 
a Q 2 = 0.0017 ± 0.0003, a\ = 0.00013 ± 0.003 . 

These were the values fitted in Ref. [fif . 



(41) 
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C. Fitting the tt-tt scattering lengths 



In Ref. ChPT was used to fit the scattering lengths: a best fit being obtained in that 
parametrisation with: 

ai = -0.0092 a 2 = 0.0080 (42) 

which yields: 

ag = 0.152, a 2 = -0.0451, a\ = 0.0363, ajj = 0.00142, ajj = -0.00109 . (43) 

Before proceeding to a calculation of the scattering lengths in the GCM it is of interest to 
determine whether a fit is possible at all and how good it might be; i.e., if there are values of 
the coefficients K\, K 2 and K% in the GCM which provide a good fit to the scattering lengths. 
This is not a well posed problem since in general, as we have seen, the scattering lengths are not 
well known. For the purpose of illustration we will then simply choose to attempt to fit a "model 
experimental data set": o[] and Og from Eq. fl40]); a\ and from Eq. (]4l"|); and a 2 from Eq. ( f43| ) 
(based on the plots in Ref. Q): 

ag = 0.20, a 2 = -0.037, a\ = 0.038, at] = 0.0017, ajj = -0.0011 . (44) 

We find a best fit (defined as that set of Ks which minimise the sum of the squares of the 
relative difference between the model experimental data and the calculated values) with 

K x = -0.000455, K 2 = 0.000414, K 3 = 0.00150 , (45) 

which corresponds to ol\ = —0.00972 and a 2 = 0.00901, and these values yield 

a[j = 0.15, a 2 = -0.042, a\ = 0.035, a° 2 = 0.0017, a\ = -0.0011 . (46) 

This fit has a mean deviation in the magnitude of the relative error of 10%. 

It should be remarked at this point that, as far as the at-threshold parameters are concerned, 
the only contribution that pion loops will make is to modify the values of the i^-parameters 
calculated from the expressions in the appendix; i.e., pion loops will only modify the expressions 
defined implicitly via Eq. ( |39| ) by requiring that 

^tree _^ ^-loops 

and hence these equations are actually form invariant and may be looked upon as providing a very 
general parametrisation of these at-threshold parameters. 

If we arbitrarily set K\ = then the expressions in Eqs. (??-|37|) reduce to exactly those of 
Ref. Q (with a minor correction of a typographical error in Tg therein). In this case the "best fit 
values" of a\ and a 2 are: 
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ai = -0.00955 a 2 = 0.00897 (48) 

which differ little from Eq. fl4^ ) and correspond to 

K t = 0.0, K 2 = 0.000844, K 3 = 0.00149 . (49) 

[Obviously, K 2 has changed in order to compensate as well as possible for the loss of the K\ degree 
of freedom: In Eq. @ K 2 - K x = 0.000869 to be compared with K 2 in Eq. (||).] Using Eq. @) 
one obtains: 

a° Q = 0.15, al = -0.045, a\ = 0.037, a° 2 = 0.0018, a\ = -0.0011 . (50) 

This fit has a mean deviation in the magnitude of the relative error of 11%. In this illustrative 
example the extra freedom associated with the coefficient K\ allows for a marginally improved fit 
in the GCM. 

We have summarised these results in Table. 1. 



IV. MODEL EVALUATION OF THE SCATTERING LENGTHS 



A. Point Meson Limit 



We turn now to the direct evaluation of the scattering lengths in the GCM using the formulae 
in the appendix that relate the coefficients K±, K 2 and K3 to the quark-gluon substructure in the 
GCM. The first and simplest calculation is the point meson limit. Writing the fermion propagator 
in the form 



G(p) = K(p 2 ) (-i 7 • P A(p 2 ) + B{p 2 )) (51) 

with 

K( V 2 ) = p 2 A ^2 + £( p 2)2 ( 52 ) 

the point meson limit is obtained with 

A(p 2 ) = 1 and B(p 2 ) = const. (53) 



This effectively reproduces the results one would obtain in the NJL [1C]. In this case f n and (qq) 
are both given by divergent momentum integrals which must be regularised but the Ks are finite 
and, in fact, 

K 1 = K 2 = K 3 = — Ut . (54) 
6 96?r 2 v ; 



(This corrects a typographical error in Ref. 0.) In this limit r n is also finite and from Eq. (91) 
one obtains: 



11 



which agrees with that in the infinite cutoff limit of the NJL model, [Ej] as it should. 

In order to present results in this case we simply take the experimental values of f n = 0.093 GeV 
and nix = 0.1385 GeV (these values can always be arranged in the NJL) and we find that: 

,V = 0.58 fin, 

ag = 0.17, al = -0.048, a\ = 0.036, a% = 0.0020, a\ = 0.0 . 

One sees that even in the (properly regularised) point meson limit the GCM provides a reasonably 
good description of the scattering lengths. 



B. Finite Size Pion 

It is simple matter to go beyond the NJL limit of the GCM. It will be recalled that in the GCM 
the fermion propagator is obtained as a solution of the fermion SDE: Eq. (|i"2|); which, in general, 
is a pair of coupled, nonlinear, integral equations the solution of which is completely determined 
once the model gluon propagator is specified. The solution of these equations in phenomenological 
models of QCD is an important area of research in its own right and has received a good deal of 



attention [14,^8],[2£,3C,31]. The characteristic features of the fermion propagator obtained in all of 
these these studies are the same and may be distilled: In the spacelike region, A is a slowly varying 
function of p 2 with A(p 2 = 0) > 1 and A(p 2 = oo) = 1; B provides most of the structure in the 
propagator with B(p 2 = 0) ~ 1 GeV and with the asymptotic form expected from the Operator 
Product Expansion |^] and QCD Renormalisation Group: 

K*)Uoo - (57) 

with A = 12/(33 — 2iVf) [we use iVf = 3 throughout], k = (In // 2 )~ A < qq > /J , a dimensionless 
renormalisation point invariant, and x = p 2 /Aq CD . 

Taking these considerations into account, and following Ref. [33| in which fits to various SDE 
solutions were made, the following model forms of A and B, that incorporate all of the physical 
information and QCD input available from the SDE studies, suggest themselves: 



MP 2 ) = ^ (58) 



n , 2^ , (, , r, ft 2\1\ P 4 47T 2 A (^[VA^cd]) A (^)lGeV 2 

B(p 2 ) = h (l - tanh [6a (& 3 + p 2 )}) + J^ — - J ' y-x ■ (59) 

p 2 In 64 + 



A 2 

QCD , 



(We remark that these fitting forms should only be used at spacelike momenta: p 2 > 0; i.e, they 
are useful for interpolating the SDE solutions but not for extrapolating them into the complex 
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p 2 plane. Studies of the solution of model fermion SDEs in the complex plane can be found in 



Refs. ||.) 

The experimentally determined QCD mass scale, Aqcd ; introduced into QCD by renormal- 
isation, can be inferred from a number of experiments and the average value for Nf = 4 in the 



modified minimal subtraction renormalisation scheme is: [35 

Aqcd = 238 ± 30 ± 68 MeV . 

For Nf = 3 this corresponds to 

Aqcd 290 MeV 

with a charm quark mass of 1.5 GeV and this is the value we use herein. 
A reasonable choice for the value of the quark condensate is 



(60) 



(61) 



PlGeV 2 



(uu) 



lGeV 2 



(dd) 



lGeV 2 



(qq) 



lGeV 2 



(0.255GeV) ; 



(62) 



which follows from Refs. [29,36| and m n is then determined by Eq. (|29| ) once f n is calculated from 
Eq. (0). [We fix (m« 2=1GeV2 +mJ 2=1GeV2 ) = 10 MeV.] 

The undetermined parameters in the SDE based model fermion propagator are 61-64. Choosing 
values for these parameters is completely equivalent to choosing a model form for the IR (small k 2 ) 
behaviour of the model gluon propagator in the GCM. This is because the UV (large p 2 ) behaviour 
of A and B is guaranteed by the requirement that the model gluon propagator manifest asymptotic 
freedom at large k 2 ; i.e., 

Xtt 

k 2 large 



F T (A; 2 



In 



fc 2 
A 2 



(63) 



This is established by the SDE studies of Refs. ||,||,||,|gj3| 
In the fit of Ref. §| 

61 = 0.7336 GeV 6 2 = 4.779 GeV~ 2 fo 



-0.1435 GeV 2 



(64) 



which were chosen to give the best possible fit to the SDE results of Ref. |37[] on the domain 
s = p 2 £ [0, 1] GeV 2 . It is worth remarking that the parameter 64 in Eq. ( |59] ) is important only in 
that it should be greater that one so that the small p 2 behaviour is governed by the "tanh" piece 
of the model. As a consequence we set 



64 



(65) 



and consider it fixed hereafter since our results are insensitive to increasing it by two orders of 
magnitude. This is also the reason for the factor p 4 /(l + p A ) in Eq. (|59|). 
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The model quark propagator obtained from Eqs. ( |58D and fl59|) can be used to calculate f n , 
m n , r n and the partial wave amplitudes. As a first step we calculated these quantities using the 



parameters of Eq. ( |64| ) which tests the model SDE solution of Ref. [37|. This yields 

/ w = 0.074GeV, m n = 0.17GeV, rv = 0.72fm, , . 

(66) 

ag = 0.37, a 2 , = -0.15, a{ = 0.072, a° 2 = -0.0030, a\ = -0.0067 . 



The value of f w agrees with that in Ref. |J7j while the other quantities have not previously been 
calculated in this model. 

These results are not generally satisfactory, a fact that is not surprising since the parametri- 
sation of the gluon propagator in Ref. |J7j was simply a "first guess"; i.e., no systematic attempt 
was made to improve the fit to static properties of mesons. The obvious next step is to determine 
whether better agreement with experiment can be obtained simply by varying bi, b 2 and 63. 

We allowed 61, b 2 and 63 to vary in order to perform a least squares fit to f % = 0.093 GeV, 
m n = 0.1385 GeV, r w = 0.66 fm, ag - a| from our model experimental data set, Eq. (Q), and 
values of Tq(x), Tq(x) and T/(x) at x = 0.28 GeV 2 taken from the plots of Ref. pL This procedure 
yields 

61 = 0.5901 GeV, b 2 = 3.0508 GeV" 2 , 63 = -0.2961GeV 2 . (67) 

The mass function obtained with this parameter set is plotted in Fig. ^ and, for comparison, 
we also plot B(x = p 2 ) obtained with the parameters of Eq. (|64|) in this figure. It will be observed 
that these functions differ only for |p|/Aqcd ^4; i.e., for infrared momenta. It is important to 
recall that Eq. fl64j) yields an exact fit to the solution of the fermion SDE obtained when the 
model gluon propagator has an integrable IR singularity oc S 4 (k). This can be interpreted p9| as 
a regularisation of the 1/k 4 behaviour argued for in Refs. [38|,39|40],[ll|]. 



The parameters of Eq. ( p7|) lead to significantly improved agreement with the data: 

f n = 0.091 GeV, m n = 0.141 GeV, r ff = 0.59fin, , . 

(68) 

ag = 0.17, a 2 = -0.052, a\ = 0.032, a Q 2 = 0.0011, a\ = -0.00085 . 

(The values of the ET-coefficients are listed in Table. I). It is of interest to note that is too small 



by ~ 10% which, consistent with the lattice analysis of Ref. [20|, allows for a small correction from 
the 7r-loop at the physical value of m n . 

It is clear now that, apart from m n , the characteristic quantities of the pionic sector are a 
measure of the IR structure of the model gluon propagator; i.e., the effective interaction between 
quarks at small k 2 . It is worth emphasising here that this follows because there is a one to one 
correspondence between the IR behaviour of the quark propagator and that of the model gluon 
propagator. In fact, in principle it is possible to invert the fermion SDE, Eq. (|l~2^), to obtain 



the gluon propagator that corresponds to Eq. (68). For our purposes herein, however, this is not 



necessary. (Although it should be remarked that the form of A and B is consistent with the 
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presence of an integrable singularity in the model gluon propagator at k 2 = 0. flll| , f27| , pS| 1 ) The one 
thing it is necessary to recall in this connection is that once one has a form for B(p 2 ) then the form 
of the pion's Bethe-Salpeter amplitude, T, follows immediately because of the identity between B 
and r in vector exchange theories with DCSB: [12] 



T^p^P 2 = 0)^B(p 2 ) . (69) 

To further illustrate the importance of the IR part of the fermion propagator we have calculated 
the characteristic parameters of the 7r-sector using Eq. ( p7| ) but suppressing the UV tail contribution 
in Eq. (|59|); i.e., neglecting the In term, in which case: 



U = 0.090GeV, r 7r = 0.58fm, , . 

' (70) 
ag = 0.18, a 2 Q = -0.055, a\ = 0.034, a£ = 0.0012, a\ = -0.00094 . 

Clearly, the UV tail only contributes a small amount. In this connection we have also studied the 
sensitivity of our results to variations in Aqcd in our parametrisation. Changes of ± 50% produce 
insignificant (< 1%) changes in the characteristic parameters. This is somewhat artificial, however, 
since in reality the coefficients 61-63 depend on Aqcd through the solution of the fermion SDE. 



We present the partial wave amplitudes of Eqs. (??-|37|), calculated with the parameter set of 
Eq. (|67D, in Fig. The agreement with the experimental data, even away from threshold (out to 
E ~ 0.7 GeV), is quite good. 

Clearly, the GCM is able to provide a good qualitative and quantitative description of the 
pionic sector of QCD with only three parameters. 



V. SUMMARY AND CONCLUSIONS 

We have studied the pionic sector of QCD using a model field theory which incorporates 
asymptotic freedom, dynamical chiral symmetry breaking (DCSB) and quark confinement: the 
Global Colour-symmetry Model (GCM). Since the GCM manifests DCSB, the results of chiral 
perturbation theory (ChPT) are necessarily contained within it. In this connection, the GCM allows 
for a microscopic interpretation of the parameters in ChPT; i.e., it relates them to properties of the 
model dressed quark and gluon propagators in the GCM which are the fundamental ingredients of 
this model. 

In deriving the 0(E 4 ) action for tt in the GCM we demonstrated explicitly that ChPT phe- 
nomenology is contained within the GCM since our action contained that of Ref. M plus an 
additional term allowed by chiral symmetry. We showed that, at the simplest level, this extra 
term, characterised by a single parameter, enables a marginally improved fit to the tt-tt scattering 
lengths. 

Calculations in the GCM are tied to the large body of work that has been completed in the 
area of phenomenological coupled Schwinger-Dyson^->Bethe-Salpeter equation (SDBSE) studies in 
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QCD and, in fact, the application of the GCM can be interpreted solely within that framework. 
In the calculation reported herein we used the SDBSE studies to motivate a form for the quark 
propagator - a form that incorporated the physical input of asymptotic freedom and the operator 
product expansion; i.e., the large p 2 (UV) behaviour of the propagator was completely determined 
- so that our model propagator had only three parameters. These parameters, which determine 
the small p 2 behaviour of the propagator, reflect the fact that little is known about the form of 
the quark-quark interaction at small k 2 (IR). (The SDBSE studies indicate that there is a one to 
one correspondence between the IR behaviour of the gluon propagator and the IR behaviour of the 
quark propagator.) 

Using the three parameter form of the GCM we calculated the quantities that characterise the 
pionic sector of QCD: f n , m n , r n and the lowest five partial wave amplitudes for ir-ir scattering at 



tree level. We allowed the parameters to vary in order to obtain the best fit possible. In Sec. IV B 
we demonstrated that the GCM provides a very good description of this sector of QCD with a 
good representation of the partial wave amplitudes well beyond threshold. 

It is worth reiterating here that in our study there are no infinite regularisation parameters, as 
there are in ChPT, even if we proceed to incorporate pion loop diagrams. We have a model, which 
can also be understood as a relativistic potential model, in which the parameters are tied to the 
behaviour of the quark-quark interaction in the infrared. 

The results we have presented herein provide evidence that a "two-point model" of QCD; i.e., 
a model in which coloured quark currents interact via dressed gluon exchange (where the dressing 
models the effect of gluon and ghost polarisation diagrams on the gluon propagator but where 
interactions that depend on explicit three- and four-gluon vertices are neglected) provides a good 
qualitative and quantitative understanding of low energy phenomena in QCD. It is worth remarking 
in this connection that a better experimental determination of the scattering lengths would allow 
for more stringent constraints to be placed upon the nature of the effective interaction between 
quarks at small k 2 . 

1 ome of the calculations described herein were carried out using a grant of computer time and 
the resources of the National Energy Research Supercomputer Center. The work of CDR was 
supported by the US Department of Energy, Nuclear Physics Division, under contract number 
W-31-109-ENG-38, with additional funding from the US National Science Foundation, Division of 
International Programs, under grant number INT-92 15223. The work of RTC was supported in 
part by the Australian Department of Industry, Technology and Commerce. 



COEFFICIENTS IN Q[U] 

In this appendix we present the expression for r n and expressions for the coefficients Ki, K% 
and K% in Eq. (??). The manner in which they are calculated is described in Ref. j/J and the 
expressions presented here correct some typographical errors therein. Further, they extend those 
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expressions by incorporating the dependence on the derivatives of A(p 2 ). 
We have 

Ki = -b + c - 2h + / + g 

K 2 = a K2 - d + 2f + g + i - j - 21 + 2m + 2n K2 
K ?> = a K3 + e-2g -2i + Ak + 2l - 2m + 2n Kz 

where 

aK * = !wf rAbA ( aA + 2 ^ AA ' + «W) + h 4{AA ' + q2[A ' ]2)2 ) ' 

aK3 = Iwf k ' bA ( aA + 2 ^ AA ' + «W) - \^ AA ' + ?W) 2 ) . 

b= l f Jl* B 2 + A 2 KK f + 

2 J (2vr) 4 V 

|-(if 2 [A'] 2 + 2AKA'K' - A 2 [K'} 2 + 2AK 2 A" + A 2 KK") + 

q\K[A'] 2 K' - AA'[K'} 2 + -K 2 A'A" + AifA'if" + -AK 2 A"') + 

2 6 



c = 



g 6 (-[A'] 2 [^'] 2 + K[A'] 2 K" + l/^A'A'")) , 

rJik B2 ( K " +q2K, " + \ qAKn 



d= \IWf {K 2 (-(BB'f + q 2 BB'(BB" - [B'} 2 )) 

+ i q A {A{KK" - K ,2 ){BB') 2 + K 2 {2B 2 B'B'" - (BB") 2 - [B'} 4 ))^j 

e = I J ^K(2K'(BB') 2 + KB 2 B'B" + 0[B'] 3 ) 

+- (/ 2 (4(ifif" - K' 2 ){BB') 2 + K 2 {2B 2 B'B"' - {BB") 2 - [Bf))) 
6 / 

i= \Sw? k2 ( (bb ' )2 + 9 2BB '([ B 'i 2 + BB "> + 59" (i B 'i 2 + BB "> 2 ) 
*=hlw^ ,K2 B f + BB ") 2 

= I / g 2 K 2 BB' ( [Bf + + -q 2 (W'B" + BB'")") 

4 7 (2vr) 4 V 3 / 

y q 2 K 2 B 2 ^AB'(A B K' +2B A' K + AB'K)+^q 2 



A ABA' B K 

+8 B [A'} 2 B' K + 5AA' [B'f K + 2AB B' A" K + AB A' B" K 
+q 2 A'(4B A' B' K' +5 A' [B'} 2 K + AB B' A" K + B A' B" K 



j 1 • \-BB'K 



(2ny 



-A K (2A'BB' + A{[B'} 2 + BB"}) + -q 2 IaABA'B'K' 
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n K 2 



n Ka 



+5B[A'] 2 B'K + 5AA'[B'] 2 K + 2ABB' A" K + ABA'B"K 

+q 2 A'(4BA'B'K' + 5A'[B'] 2 K + 4BB'A"K + BA'B"K)^ ) (84) 

fc = l/(0 ?W[B ' ]3 (85) 

l = \SWf {^B 2 [B']\[B'\ 2 + BB") + i i^i? 3 [£?) (86) 

m = \lwf ^ K ' B 'l B T) ( 87 ) 

/ (^ tf* 4 * W ( A2 + f ?W + (88) 

/ (^ ( A2 ~ I ^ AA ' + (89) 

In these expressions 

F> = jL F {q 2 ) (90) 

with F" representing the second derivative, etc. 

The 7r charge radius obtained from the generalised impulse approximation diagram, Fig. |2[ is 



W (2t) 4 
where 

H(x = q 2 ) = B (-A 2 B 2 B' K" K 2 x 2 - A^B" K 1 K 2 x 3 - 36 A 2 B 3 K' K 2 - A^B' K" K 2 x 3 
+A 2 B A 1 A" K 3 x 3 - 30 A B 3 A' K 3 - 42 A 2 B 2 B' K 3 - 24 A 2 B 3 K" K 2 x 
-48 A 4 B K' K 2 x - 54 A B 3 A' K' K 2 x - 42 A 2 B 2 B' K' K 2 x - 72 A 3 B A' K 3 x 
-18 B 3 A' 2 K 3 x - 24 A*B' K 3 x - 42 A B 2 A' B' K 3 x - 24 A B 3 A" K 3 x 
-21 A 2 B 2 B" K 3 x - 4 A 2 B 3 K' 3 x 2 + 4 A 2 B 3 K" K' K x 2 - 6 A B 3 A' K' 2 K x 2 
-2 A 2 B 2 B' K' 2 K x 2 - 24 A 4 B K" K 2 x 2 - 3 A B 3 A' K" K 2 x 2 - 68 A 3 B A' K' K 2 x 2 
-2 B 3 A' 2 K' K 2 x 2 - 28 A 4 B' K' K 2 x 2 - 18 A B 2 A' B' K' K 2 x 2 -8 A B 3 A" K' K 2 x 2 
+3 A 2 B 2 B" K' K 2 x 2 -20A 2 B A' 2 K 3 x 2 - 40 A 3 A' B' K 3 x 2 -UAB A' B' 2 K 3 x 2 
+2 A 2 B' 3 K 3 x 2 - 36 A 3 B A" K 3 x 2 - 5 B 3 A' A" K 3 x 2 -4A B 2 B' A" K 3 x 2 
-8 A A B" K 3 x 2 + 2 A B 2 A' B" K 3 x 2 + 4A 2 B B' B" K 3 x 2 - 4 A*B K' 3 x 3 
+4 A^B K" K' Kx 3 - 6 A 3 B A 1 K' 2 K x 3 -2 A 4 B' K' 2 K x 3 - 3 A 3 B A' K" K 2 x 3 
-WA 2 B A' 2 K' K 2 x 3 - 10 A 3 A' B' K' K 2 x 3 - 4 A 3 B A" K' K 2 x 3 -6AB A >3 K 3 x 3 

(2 /I' 2 D / t^3„3 o ^3 t}' A>l TsSJA o /i3 aI u" tz^JS 



-6 A 2 A"B' K 3 x 3 - 2 A B' A" K 3 x A - 2 A 3 A' B" K 3 x 3 .j (92) 
which is obtained with the fermion-photon vertex Ansatz of Ref. 42]: 
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V^q) = \[A{p)+A{q)\^ (93) 
+ (P + 



p2 _ g2 



|[A(p) - A( g )]i( 7 ■ p + 7 • ?) - *[5(p) - 5(g)] 



which ensures local current conservation in this calculation. 

Computer code for all of the quantities listed in this appendix in either MATHEMATICA or 
FORTRAN format is available upon request. 
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FIGURES 

FIG. 1. This quark loop diagram contributes to f w . The filled circles represent {Tr\qq) 
Bethe-Salpeter amplitudes, the thick external lines represent the tt field and the thin internal 
lines represent dressed quark propagators. 

FIG. 2. This diagram provides the piece of r n , Eq. (|l|), that is regular in the chiral limit. The 
thick, straight external lines represent the incoming and outgoing tt, the filled circles at the tt legs 
represent the (ir\qq) Bethe-Salpeter amplitudes, the wiggly line represents the photon, 7, the shaded 
circle at the 7 leg represents the regular part of the dressed quark-photon vertex (which satisfies 
the Ward Identity), Eq. (|93"1), and the thin internal lines represent the dressed quark propagator in 
the model. 



FIG. 3. The function B(p 2 ) of Eq. ( |59| ) obtained with the parameters in Eq. (|67|): solid line; 
compared with that obtained with Eq. (|64|): dashed line. 



FIG. 4. The partial wave amplitudes of Eqs. (??-|37|) obtained with the best fit described in con- 
nection with Eq. (j38|): solid line. In the plots of T®(x), If(x) and Tq (x) we show, for comparison, 



the results that would be obtained with the Weinberg result for A(s, t, u) in Ref. |22j (short-dashed 
line): this form gives T 2 °(x) = = If (jc). In these plots we also show the extrapolation to threshold 
of Roy equation fits to high-energy tt-tt scattering data for these partial waves: long-dashed line. 
The data are extracted from Fig. 4 in Ref. . We have defined x = E 2 — Am 2 . 
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TABLES 

TABLE I. In this table we summarise our results for the scattering lengths calculated in 



Sees. Ill B and IV. The column labelled "Exp" lists the "model experimental data" set that 



we constructed; see Eq. (44). Fit 1 is the best fit to "Exp" possible (in the sense of minimising the 
sum of the squares of the differences) in the GCM when the Ks are allowed to vary; see Eq. (|46|) . 
Fit 2 is the best fit possible when one sets K\ = and allows 1^2 3 to vary and corresponds to the 
best fit to "Exp" possible with the parametrisation of Ref. [|j; see Eq. (|50|). Cal. 1 is the calculated 
result obtained in the point meson limit of the GCM and corresponds to the result that would 
be obtained in the Nambu-Jona-Lasinio model [l/(967r 2 ) = 0.00106]; see Eq. fl56|). Cal. 2 is the 
calculated result using the quark propagator given by Eqs. ( [58| ) and (|59|) and using the parameters 
of Eq. (|3); see Eq. (I 
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